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The inverse scattering theory is a basic tool to solve linear differential equations and some Partial 
Differential Equations (PDEs). Using this theory the Korteweg-de Vries (KdV), the family of 
evolutionary Non Linear Schrodinger (NLS) equations, Kadomtzev-Petviashvili and many more 
completely integrable PDEs of mathematical physics are solved, using Zacharv-Shabath scheme. 
This last approach includes the use of a Lax pair, and has an advantage to be applied to wider 
class of equations, like difference equations, but has a disadvantage to be used only for "rapidly 
decreasing solutions". This technique is also intimately related to completely integrable systems. 
The identifying process of a Lax pair, a system and finally the "scattering data" is usually a difficult 
process, simplified in many cases by physicals models providing clues of how the scattering data 
should be chosen. In this work we show that the scattering data can be encoded into singularities 
of a very special mathematical object: J-unitary, identity at infinity matrix-valued function, which, 
if evolved once, solves the inverse scattering theory, if evolved twice solves evolutionary PDEs. The 
provided scheme seems to be universal in the sense that many (if not all) completely integrable PDEs 
arise (or should arise) in this manner by changing the so called "vessel parameters" (for examples, 
solutions of KdV and evolutionary NLS equations are presented). The results presented here allow 
to study different fiows (commuting and non-commuting) in a unified approach and provide a rich 
mathematical arsenal to study these equations. The results are easily generalized to completely 
integrable PDEs of n (< oo) variables. 

PACS numbers: 02.30.Jr, 03.65.Nk, 02.30.Ik 
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II. DESCRIPTION OF THE METHOD 
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A standard technique to solve some non linear PDEs, 
initiated in (CG67j and gen eralized to a Zacharov- 
Shabbath scheme in [ZABS74| includes for example the 
Korteweg-de Vries equation 
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(1) 



where q^, denote the partial derivatives. This equation 
arises in the study of waves on shallow water surfaces. A 
second important example is the Non Linear Schrodinger 
(NLS) evolutionary equation: 



»/3t+/3xx + 2|/3p/3 = 0, (3{x,0) = (3{x), 



(2) 



where /3 = (3{x,t) is a complex valued function of two real 
variables x,t and P(x,0) = I3{x) is the initial condition, 
defined usually on R. This equation plays a special role in 
optics and water waves, we show how these two impor- 
tant equations of mathematical physics arise as special 
cases of a theory of vessels, recently develo ped by the au- 
thor, initiated by works of M. Livsic |Ls01 [. As we would 



like to focus on the technique and not on mathemati- 
cal aspects of the theory, we will not actually define the 
vessel itself, but will use all its necessary ingredients to 
construct solutions of many completely integrable PDEs. 



A. Scattering theory 

Suppose that p e N and let 

cri = crj, 0-2 = cr2, 7 = -7*, 

be three p x p matrices such that ai is invertible. We 
call such a triple as vessel parameters. For example, 
some triples for the case p — 2 are given in Table H] By 
a scattering data we mean a p x p function 5(A), called 
scattering data matrix possessing a realization 

S{X) =1- B*Xo i(A/ - A)-^Baai, 

AXo + XqA* + BocriBo = 0, X*^ ^ Xq (3) 

where for an auxiliary Hilbert space Ji the operators act 
as follows: Xo,A : n H, Bp - CP H. Such re- 
alizations where studied in jBGRQOl l and we present the 



TABLE 1. Parameters corresponding to Sturm-Liouville (SL), 
Non Linear Schrodinger (NLS) Equation and Canonical sys- 
tems. / - the identity matrix, - the zero matrix. 
P=2 
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regular case, where all the operators are bounded to clar- 
ify explanations and to simplify formulas, but the general 
case with an unbounded A is presented later in Section 
mil Letting the operators evolve with x as follows: 



dx 



B{x)ai = -AB{x)a2 - B{x)-f, B(0) ^ Bq, (4) 
[x) = B{x)(j2B*{x), X(0) = Xo (5) 



d . 
dx 



we obtain that for the self-adjoint X(x) the Lyapunov 
equation 

AX{x)+X{x)A* +B{x)aiB{x)=0 (6) 
holds (see jAMll) and the function 

^(A.a;) = I - B*{x)X-\x){XI - A)-^B{x)ai (7) 
serves as a Backlund transformation: 



Theorem II. 1 r |Mel09LlA]V[12l] ^ Fix A ^ spec (A) 
and let u{X, x) be a solution of the input Linear Differ- 
ential Equation (LDE) 

d 

- ai—u{X,x) + {a2X + l)uiX,x) =0. (8) 

Then the function y{X,x) = S{X,x)u{X,x) is differen- 
tiate and satisfies the output LDE 

d 

- + ('^sA + -i^{x))y{X,x) = 0, (9) 

where using Hq{x) — B* {x)%r^{x)B(x) 

7*(2^) = 7 + (y2Ho{x)ai - aiHo{x)a2. (10) 

In order to explain this Theorem, the Table U is useful. 
Taking the vessel parameters cti , (72 , 7 from a row in this 
Table and using an arbitrary function 5'(A) as in ([3]), we 
will obtain that S'(A, x), defined in ([7]) maps solutions 
of a simple LDE, with constant coefficients, to a more 
complicated one, defined by 7*(a;). This is an analogue 
of the scattering operator. Moreover, in SL case it is 
also possible to characterize precisely the equation 

2/1 (A, a;) 
2/2 (A, a;) 

second entry y2{X,x) is uniquely derived from yi{X,x), 
and the first one satisfies the Sturm-Liouville LDE 



It turns out [Mel09l| that for y(A, 



the 



dx^ 



2/1 + qix)yi{x) = -iXyi, 



where g(x) = -2^ ln(dct X^ ^X(a;)) (see |Mel09l | for de- 
tails). Worth noticing that the famous tau- function is 
defined in general on the basis of this example as follows: 



B. Evolutionary equations of SL-KdV type 

An interesting phenomenon occurs when we let the op- 
erators evolve with respect to a variable t. We choose to 
represnt the following two evolutions: 

^^Bix,t) ^ iA^B{x,t), B{x,0)=B{x), (12) 



d_ 

dt 



X^iABa2B* -iBa2B*A* +iB-fB*, (13) 
X(a;,0) = X(a;), 



which in the case of SL parameters create solutions of the 
KdV equation ([T]). That's why we also call these evolu- 
tionary equations as of SL-KdV type. In the case when 
all the operators are bounded it is a matter of simple 
calculations to check the following 



Theorem II. 2 ([Meldll^ For the SL vessel parame- 
ters, let B{x,t),X{x,t) be two operators satisfying 
dH), ©, (dll), dm). Then the function q{x,t) = 

-2^1n(detX(7^X(x)) satisfies 

Theorem II. 3 ( (Meldl l For the NLS vessel param- 
eters, let B{x,t),X(x,t) be two operators satisfying 
(HI), ©, ([HD, 1131). Then the function P{x,t) = 

satisfies ^ for ^^{x,t), defined by 



l]i.,{x,t) 



If the reader is interested in what happens for the canon- 
ical systems parameters, the answer is as follows and its 
proof is omitted. In order to understand how it is proved, 
see Section III El 

Theorem II. 4 For the canonical systems vessel param- 
eters, let B{x,t),X{x,t) be two operators satisfying (j4]), 

-2il3{x,t) ih{x,t) 



dSD, (ini), (USD. Then -f,{x,t) 



with h{x,t) 



(x + K)^ 
valued function j3{x,t) satisfies: 



ih{x,t) 2il3{x,t) 
4/32 (K eM.) and the real 



d 



[x + Kf 



{x + Kf 



(14) 



t{x) = det(Xo iX(a;)) 



(11) 



Finally, the class of potentials, corresponding to the reg- 
ular vessels is as follows: 

Theorem II. 5 Suppose that B{x),X{x) are defined by 
([4]), ([5]) using a bounded operator A, then 7* (a;) is in- 
finitely differentiable at each point, where t{x) 7^ 0. 

Proof: Notice that using (jj]), ([5]) the expression pO)) is 
infinitely times differentiable. □ 
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C. Commuting flows. Hierarchies 



The formulas (|T2)) . ([T3| have a special property. It is al- 
most readable that the mixed second derivatives, applied 
to B{x,t) are equal: ^^B{x,t) = ^g^B{x,t). But the 
same is true for lL{x,t), which is shown in |Melbl | : 

Lemma II. 6 Suppose that B{x,t),X{x,t) satisfy 
m, (HI, (HI. Then 



92 



dxdy 



ix,t) - 



92 



dydx 



X{x,t). 



Similar formulas may be used to derive evolutionary 
PDEs, constructed from other commuting flows. Let us 
demand 



-{Bai)^AY^A'Bm, 

i=0 



rrin = (-I)"m^ 



and 



where 



i=0 



Y» = J2^-iyA''-^Bm,,B*iA*y 



(15) 



(16) 



(17) 



i=0 



It was obtained in [Melbl . Lemma 14] that in a basic case 
the corresponding flows commute: 

Theorem II. 7 Suppose that B{x,t) satisfies (|4|) and 

B[ = {^/^ArB',. (18) 
Suppose also that X(x,i) satisfies ([S]) and (I16|) . then 

dxdy dydx' dxdy dydx 

The same result conjecturally holds for dTSl) . (ITCl) . when 
the polynomial with matrix coefficients po{X) = Ao'2 + 7 

n 

commutes with ^ X^mi. In a particular case of com- 
muting flows it was also discovered jMelbl . Theorem 19] 
the corresponding family of PDEs, which also constitute 
a part of the KdV hierarchy: 

Theorem II. 8 Suppose that B{x,t),X{x,t) satisfy ([!]), 

1 3 

(0, dSD and (dH). Assume that bo = --fixxx + -j^il^xf 

and define recursively a differential polynomial in P{x) 
as follows (n = 0,1,2 .. .) 

4(&„+l)2; = -i{bn)xxx + ^i(l3xbn)x- (19) 

Then the 1,2 entry of jf{x,t), defined by (llOp satisfies 

coinciding for n = with the ( integration with respect to 
x of) the KdV equation ([T]). 



D. Completely integrable systems corresponding to 
SL-KdV type 

Let (u(x, i)/x(x, i)/y(a;, t)) be the (input/ state/ out- 
put) triple of the following system of equations 

^x(i, x) — AyL{t, x) 4- B{x, t)criu(i, x), 
^x{t,x) = B{x,t)a2u{t,x), (20) 
y{t,x) = u{t,x) - B*{x,t)X-'^{x,t)x{t,x). 

which is overdetermined, namely only u{t, x) satisfying 
is allowed. Then the output satisfies 

r 9 d / M / N f, 

[-CTi^ +cr2^ +-/*{x,t)\y(x,t) = 0, 

which can be checked using some algebraic manipula- 
tions analogously to the proof of Theorem lII.il (Backlund 
transformation). This system is integrable (for legal in- 
puts!), because the mixed second derivatives are equal: 



d^ 



-x(t,x) 



d^ 



-x{t, x) 



dxdt ' dtdx 

which can be checked by the use of system equatio ns pop , 
the equations (g]), (P and dHJ (see |MLsV9511LsOH for de- 
tails) . Performing here separation of variables u(t, x) = 
e^^u{X,x), x{t,x) = e*^x{X,x), Y{t,x) = e^^y{X,x) we 
wiU arrive to the differential equations dH), ([9]) for u{X, x), 
y{X,x) respectively. So, vessels indeed come from com- 
pletely integrable systems by separation of variables. 

E. Evolution of ^(A, x) for SL-KdV type 



Notice that the Backlund transformation Theorem lII.il 
is equivalent to the fact that S satisfies 



d_ 

dx 



S = a^\a2X + -i,)S - Sn-^a^X + 7) (21) 



Moreove r, for t he equations of SL-KdV type it is possible 
to show [Melbl . (33)] that 



(22) 



where is the zero moment, while a general moment is 
defined as i/„ = B*%-'^A''B. Finally calculating 

^7* = Tri(T2B*X-^Bai - —aiB*\-^Ba2 
at dt dt 

using (jl2p . (I13p and then rearranging the terms, we will 
obtain that 

{l*)'t = -il*{Ho)'x^i+iai{Ho)lx^i+iai{Ho)'xl* (23) 

It is a matter of a little effort (after studying the mo- 
ments) to show that (|23p is equivalent to ([T]) in the SL 
case, to ([U in for the NLS vessel parameters, and to (fH)) 
for the canonical system parameters. 
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F. Higher dimensional completely integrable PDEs 

It is a very interesting question which PDEs are solved, 
if we choose to use the fam ily of N LS equations, defined 
by A.P. Fordy, P.P. Kulish |AF83| . These also constitute 
examples of applicability of the suggested scheme to solve 
the scattering theory for different values of p. 

Another important exam ple is the Kadomtzev- 
Petviashvili equation |BK7Cll |. It is a generalization of 
the KdV equation ([1]). An example of construction 
of solutions for this equation is as follows. Define for 
B = B{x,t,y) 



Bf — B^ 



Bx, 



and B ^, — B^ , Xj^ — . Then it follows using ideas in 
[Meldl l that 



(7* 



d 



which turns to be equal to {'j*)xx, resulting in KP equa- 
tion, when all the terms are transfered to the left hand 
side. 



III. THE NON REGULAR CASE 



holding for all A ^ spec(A). The equation ([5]) remains 
the same and the Lyapunov equation ^ becomes 



X{x)R*{-X) + R{X)X{x)- 

RiX)B{x)aiB* {x)R* {-X) 



0. (25) 



Notice that (HH), (gl]) coincide with (g]), (O after cancel- 
lations when A is bounded. In this case, it is possible to 
show that S{X, x) defined in ([7]) will map solutions of ([8]) 
to (0), defined by 7*, which usually fails to be differen- 
tiable, generalizing Theorem III. 5 1 



IV. STANDARD MODELS 



a. Solitons The solitons are obtained if we consider 
so called finite dimensional realizations (dimH < 00). In 
these cases the scattering data matrix 5(A) is rational. 
For the Sturm-Liouville case see ^Mel 09] . 

b. Spectrum on a curve If we fix a symmetric with 
respect to the imaginary axis Jordan curve F and consider 



The same formulas can be used when we use an un- 
bounded operator, giving us a richer family of "general- 
ized potentials" 7,. In this case the operator A : D{A) 
H is assumed to have a dense domain D{A) C H and 
constitute a generator of a Cq semi-group. Let us denote 
its resolvent by i?(A) — {XI — A)~^, then for B{x) we 
demand 

d 

= —[R(X)B(x)]ai+AR(X)B(x)a2+R(X)B(x). (24) 
ox 

Notice that we obtain regularity assumptions, neces- 
sary for the equation (|24p to be well defined: 

R{X)B{x)a2 e D{A), R{X)B{x)j € H, 



H = L'{T) = {f{^,)\ / |/(A.)|^<oo} 



and A — ifj, - multiplication by the variable of F. Then 
there are explici t form ulas for the operators B{x) — 
B{ii,x), X(a;) in |Mela| . showing that ^(A,^) will have 
jumps on F only. 

c. Discrete spectrum Choo sing % — 1^ and A = 
diag(ifc^) one construct in [Melaj S{X,x) whose poles are 
precisely at the points fc„. 

d. Fadeyev inverse scattering Using the ideas of the 
continuous and th e discre te spectrum the Fadey ev inv erse 
scattering theory (Fad74 1 was implemented in [Melaj | . 
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